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5_i ' The plane-symmetric inhomogeneous cosmological models of perfect fluid dis- 

Oh< tribution with electro-magnetic field is obtained in the framework of Lyra's 

' geometry. To get the deterministic solution, I have considered A — f{x)v{t), 

B = g{x)iJL{t) and C = h{x)fj,{t), where A, B and C are metric coefficients. It 
has been found that the solutions generalize the solution obtained by Pradhan, 
Yadav and Singh (2007) and are consistent with the recent observations of type 
la supernovae. A detailed study of physical and kinematical properties of the 
models have been carried out. 
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. 1. Introduction and motivations 

in 

In recent years, our knowledge of cosmology has improved remarkably by var- 
. ious experimental and theoretical results. The universe is spherically symmetric 

' and the matter distribution in it is on the whole isotropic and homogeneous. 

. But during the early stage of evolution, it is unlikely that it could have such 

a smoothed out picture so I consider plane symmetry which provides an op- 
protunity for the study of inhomogegeneity.The present study deals with plane 
symmetric inhomogeneous models within the framework of Lyra's geometry in 
presence of electromagnetic field. The essential difference between between the 
' cosmological theories based on Lyra's geometry and Riemannian geometry lies 

in the fact that the constant vector displacement field /3 arises naturally from 
the concept of gauge in Lyra's geometry, where as the cosmological constant A 
was introduced in ad hoc fashion in the usual treatment. Currently the study 
of gauge function and cosmological constant have gained renewed interest due 
to their application in structure formation in the universe. 

Einstein introduced his general theory of relativity in which gravitation is 
described in terms of geometry of space time. Einstein's idea of geometriz- 
ing gravitation in the form of general theory of relativity inspired the idea of 
geometrizing other physical fields. Shortly after Einstein's general theory of rel- 
ativity Weyl [1] suggested the first so-called unified field which is a geometrized 
theory of gravitation and electromagnetism. But this theory was never taken 
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seriously because it was based on the concept of non-integrability of length 
transfer. Lyra [2] proposed a modification of Riemannian geometry by intro- 
ducing a gauge function which removes the non-integrability condition of the 
length of a vector under parallel transport. In consecutive investigations Sen 
[3, Sen and Dunn [3] proposed a new scalar-tensor theory of gravitation and 
constructed an analog of the Einstein field equations based on Lyra's geome- 
try. It is thus possible [3] to construct a geometrized theory of gravitation and 
electromagnetism much along the lines of Weyl's "unified" field theory without, 
however, the inconvenience of non-integrability length transfer. 

Halford ^ has pointed out that the constant vector displacement field 0i in 
Lyra's geometry plays the role of cosmological constant A in the normal general 
relativistic treatment. It is shown by Halford [6] that the scalar-tensor treat- 
ment based on Lyra's geometry predicts the same effects, within observational 
limits, as the Einstein's theory. Several authors Sen and Vanstone [3, Bhamra 
[S], Karade and Borikar [5], Kalyanshetti and Wagmode [TU], Reddy and Inna- 
iah [TT] , Beesham [T^] , Reddy and Venkateswarlu [13] , Soleng [H] , have studied 
cosmological models based on Lyra's manifold with a constant displacement 
field vector. However, this restriction of the displacement field to be constant 
is merely one of convenience and there is no a priori reason for it. Beesham 
[13] considered FRW models with time dependent displacement field. He has 
shown that by assuming the energy density of the universe to be equal to its 
critical value, the models have the fc = — 1 geometry. Singh and Singh [T6] — 
[19], Singh and Desikan [20] have studied Bianchi-type I, HI, Kantowaski-Sachs 
and a new class of cosmological models with time dependent displacement field 
and have made a comparative study of Robertson- Walker models with constant 
deceleration parameter in Einstein's theory with cosmological term and in the 
cosmological theory based on Lyra's geometry. Soleng [14] has pointed out that 
the cosmologies based on Lyra's manifold with constant gauge vector (j) will 
either include a creation field and be equal to Hoyle's creation field cosmology 
[H] — [13] or contain a special vacuum field which together with the gauge vector 
term may be considered as a cosmological term. In the latter case the solutions 
are equal to the general relativistic cosmologies with a cosmological term. 

The occurrence of magnetic fields on galactic scale is well-established fact 
today, and their importance for a variety of astrophysical phenomena is gener- 
ally acknowledged as pointed out by Zeldovich et al. [24]. Also Harrison [25] 
has suggested that magnetic field could have a cosmological origin. As a nat- 
ural consequences, we should include magnetic fields in the energy-momentum 
tensor of the early universe. The choice of anisotropic cosmological models in 
Einstein system of field equations leads to the cosmological models more general 
than Robertson- Walker model [55] . Strong magnetic fields can be created due to 
adiabatic compression in clusters of galaxies. Primordial asymmetry of particle 
(say electron) over antiparticle (say positron) have been well established as C 
P (charged parity) violation. Asseo and Sol _27^ speculated the large-scale inter 
galactic magnetic field and is of primordial origin at present measure 10~* G 
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and gives rise to a density of order 10~^^ gcm~^ . The present day magnitude of 
magnetic energy is very small in comparison with the estimated matter density, 
it might not have been negligible during early stage of evolution of the universe. 
FRW models are approximately valid as present day magnetic field is very small. 
The existence of a primordial magnetic field is limited to Bianchi Types I, II, 
III, VIq and VIIq as shown by Hughston and Jacobs [55] • Large-scale magnetic 
fields give rise to anisotropics in the universe. The anisotropic pressure created 
by the magnetic fields dominates the evolution of the shear anisotropy and it 
decays slower than if the pressure was isotropic [IHl ISO]- Such fields can be 
generated at the end of an inflationary epoch |31] — [55]. Anisotropic magnetic 
field models have significant contribution in the evolution of galaxies and stellar 
objects. Recently, Pradhan et al. [34], Casama et al. [35], Rahaman et al. [36] . 
Bali and Chandani [37], Kumar and Singh '3F, Singh [39] and Rao, Vinutha 
and Santhi [JO] have studied cosmological models based on Lyra's geometry in 
various contexts. Later on Rahaman et al. [HHH], Pradhan et al. [331133] ob- 
tained some inhomogeneous cosmological models in Lyra's geometry. Motivated 
by these researches, in this paper, I have studied plane-symmetric inhomoge- 
neous cosmological models in presence of magnetic field with the framework of 
Lyra's geometry and also discussed the thermodynamical behaviour of universe. 

2. The Metric and field equations 

We consider the plane-symmetric metric in the form 

ds^ =A'^{dx^ ~df) + B^dy^ + C^dz^, (1) 

where A, B and C are functions of x and t. The energy momentum tensor is 
taken as 

= ip + p)u.,v' +pgi+Ei, (2) 

where p and p are, respectively, the energy density and pressure of the cosmic 
fluid, and Ui is the fluid four-velocity vector satisfying the condition 

u'u, = -1, u'x, = 0. (3) 

In Eq. Ef is the electromagnetic field given by Lichnerowicz [45] 



hih'- ( UiU^ + ^ffj) ~ ^ih' 



(4) 



where /i is the magnetic permeability and hi the magnetic flux vector defined 

by 

h^ = ^*F,,u^, (5) 
where the dual electromagnetic field tensor *Fij is defined by Synge [46] 



'Fij - -^f^ijkiF^^- (6) 
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Here Fij is the electromagnetic field tensor and eijki is the Levi-Civita tensor 
density. 

The co-ordinates are considered to be comoving so that ^ Q ^ v?' — and 
= ^. If we consider that the current flows along the z-axis, then F12 is the 
only non- vanishing component of Fij. The Maxwell's equations 



(7) 



47rJ' 



(8) 



require that F12 is the function of x-alone. We assume that the magnetic per- 
meability is the functions of x and t both. Here the semicolon represents a 
covariant differentiation. 



The field equations, in normal gauge for Lyra's manifold, obtained by Sen 
[4] as 

13 3 

% - + 2^i^] - -^9ijM'' = -87rTy , (9) 

where 4>i is the displacement vector defined as 

0, = (0,0,0, (10) 
and other symbols have their usual meaning as in Riemannian geometry. 



For the line-element HI, the field Eq. ^ with Eqs. ^ and lead to the 
following system of equations 
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where the sub indices 1 and 4 in A, B, C and elsewhere denote ordinary differ- 
entiation with respect to x and t respectively. 

3. Solutions and field equations 
Equations (HI]) - ^ lead to 

A.4 \ B44 A4 / B4 C4 \ B4C4 
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(17) 



/2S2 B B C C ' 

Eqs. (1101) - (|14p represent a system of five equations in seven unknowns A, i?, 
C, p, p, /3 and i^i2. For the complete determination of these unknowns one 
more condition is needed. As in the case of general-relativistic cosmologies, 
the introduction of inhomogeneities into the cosmological equations produces a 
considerable increase in mathematical difficulty: non-linear partial differential 
equations must now be solved. In practice, this means that we must proceed 
either by means of approximations which render the non-linearities tractable, 
or we must introduce particular symmetries into the metric of the space-time 
in order to reduce the number of degrees of freedom which the inhomogeneities 
can exploit. In the present case, we assume that the metric is Petrov type-II 
non-degenerate. This requires that 

B\\ B44 + 1B\4 \ I Cii + C44 + 2C 
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Let us consider that 



A = f{x)u{i), 
B^g{x)fiit), 

C = h{x)n{t). (19) 



Using dUl) in dH]) and we get 



g2_ hl_ 2^4 

^TT^ = = b (constant) (20) 
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and 

g_h 



91 hi 
a h 



= 2 f ^ - ^ ) = L (constant). (21) 



f^n{gh)Kb^O (22) 
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(24) 
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Eq. (gni) leads to 
and 

^ = mv'b^, (23) 

where m and n are constants of integration. 
From Eqs. ^ and we have 

1 511 /^i±^\ _ 2 /g| (2 + b) gifei 

6 <7 6 y /i b\g^ h^J b gh 

and 

2 / ^44 ^ ^4' 

Let us assume 

g = e^+^, ;. = e^-^. (26) 

Eqs. (gni) and (gS]) lead to 

W'l = M exp (^La; + [/^ ^ (27) 

where Af is an integrating constant. From Eqs. (|23p . p5|) and p6|) . we have 

2 + b\ 4 , / 2(2-6) 

b J b \ 

ML exp {lx + + 2M2 exp ^2Lx + ^^^^-^^^^ = «■ (28) 

Eq. dU]) leads to 

^-2(53^' 
Eq. (HSl) and eq. ^ lead to 

W ^ Mx + \ogN, (30) 

where iV is the constant of integration. 
Eq. (HH) leads to 



t cosh2 (^1 a \t + to) when afc < 
£cos^(y^t + to) whena6>0 (31) 

(cit + to)^ when a6 = 

where a = a6, t is constant and ci, are constants of integration. Now we 
consider the following three cases. 
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4. Case(i): ab < 

In this case we obtain 



/ = nexp(^), (32) 
^ = £cosh5(y[77|t + to), (33) 



= rcoshT^(VTa|i + io), (34) 
g^Ncxp(-^^+Mx), (35) 



2(6-2) 
1 / L6a: , 
'^^]V'"H^(^"^^'''' ^^^^ 



where r = (— ) ' 
Therefore, we have 



A = Eo exp ( J^^^^ cosh''^'" (\/Mt + to), (37) 



Lbx \ 1 

B = Goexp(^jj—^+Mx\cosh^{y^\t + to), (38) 



C = i/oexp { 2{b-2) ~ cosh^^^i + ^o), (39) 
where Eq = nr, Gq = N£, Hq = j^. 

After using suitable transformation of coordinates, the metric ([1]) reduces to the 
form 



ds^ = Eq exp 



((Sy) -osh^{VW\T){dX' - dT') + 



exp (^-(^^ + 2M^ j cosh{^/\ir\T)dY^ + 

exp (^^^-^ - 2MX^ cosh{^/\~^\T)dZ^ , (40) 

For the specification of displacement vector f3(t) within the framework of Lyra 
geometry and for realistic models of physical importance, we consider the fol- 
lowing two cases by taking /3 as constant and also /3 as function of time. 

4.1. When (5 is a constant i.e. (5 = (5^ (constant) 

Using Eqs. ([571) . ([55]) and (P^)) in Eqs. pUj) and ((T^ the expressions for pressure 
p and density p for the model (j40p are given by 

1 f2LX\ 2^ / ^ — \ 

87rp= exp l^^—^j cosh " [y\a\T'jx 
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From Eq. (|17l) the non-vanishing component F12 of the electromagnetic field 
tensor is obtained as 
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From above equation it is observed that the electromagnetic field tensor in- 
creases with time. 

The reality conditions (Elhs [47] ) 

(i)p + p>0, {ii)p + 3p> 0, 

lead to 
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The dominant energy conditions (Hawking and Ellis jJS]) 

ii)p-p>0, (m)p + p>0, 

lead to 
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The conditions and (|46p impose the restriction on /Sq. 
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4.2. When (i is a function of t 

In this case to find the exphcit value of displacement field we assume that 
the fluid obeys an equation of state of the form 



P = IP, 



(48) 



where 7(0 < 7 < 1) is a constant. Here we consider three cases of physical 
interest. 

4.2.1 Empty universe 

Let us consider 7 = 0. In this case p — p — 0. Thus, from Eqs. ([TT|) and (|14l) 
we obtain 



f2LX\ 



cosh 



■ tanh^ ^■\/| a \T ) — | a 



Halford ^ has pointed out that the constant vector displacement field in 
Lyra's geometry plays the role of cosmological constant A in the normal general 
relativistic treatment. From Eq. (HHl) . it is observed that the displacement 
vector j3{t) is a decreasing function of time. 

4.2.2. Zeldovich universe 
Let us consider 7 = L In this case 



(50) 



Using Eqs. ([38]) and ^ in Eqs. (dU and dH]), we obtain 
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From Eq. (jSip , it is observed that displacement vector /3 is decreasing function 
of time. The expressions for pressure p and energy density p are given by 
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The reality condition (EUis 1973) 

(i)/9 + p>0, (m)/9 + 3p>0, 
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4.2.3. Radiating universe 



Let us consider 7 = i In this case 
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In this case using Eqs. ([BT)) - ([39l) in Eqs. (jlip - (|14|) . the expressions for 
p and p are obtained as 
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From Eq. ()55p . it is observed that displacement vector /3 is decreasing function 
of time. The reahty conditions (Elhs [17]) 

(i)/9 + p>0, (m)/9 + 3p>0, 

and the dominant energy conditions (Hawking and Ehis |48j ) 

(j)p-p>0, (m)p + p>0, 

lead to 
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respectively. 

The expressions for the expansion 0, shear scalar cr^, deceleration parameter q 
and proper volume V for the model PO]) are given by 
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Figure 1: The plot of deceleration parameter (q) vs. time (T). 
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From Eqs. ([5^ and ([SO)) we obtain 
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The Hubble parameter iJ is given by 

LX 
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The sign of q indicates whether the model inflates or not. A positive sign 
of q corressponds to standard decelerating model where as the negative sign 
— 1 < g < 0, indicates inflations. Recent observation shows that the deceleration 
parameter of universe is in the range —1 < q < and the present universe is 
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undergoing an accelerated expansion [321 [SO] . This behaviour is clearly shown in 
Fig. 1, as a representative case with appropriate choice of constants and other 
physical parameters using reasonably well known situation. Also the current 
observation of SNe la and CMBR favour an accerating model (g < 0). From 
Eq. (pT|) . it can be seen that the deceleration parameter q < when 



T < — , tanh ^ i / — 

V (2 - 

It follows that our models of the universe are consistent with recent observations. 
Generally model (40) represents expanding, shearing and anisotropic universe 
in which the flow vector is geodetic. 

4.3. Thermodynamical behaviour and entropy of universe 

From the thermodynamics |511 152] , we apply the combination of first and second 
law of thermodynamics to the system with volume V. As we know that 

TdS ^ d{pV) + pdV (65) 

where T, S represents the tempreture and entropy respectively. 
Eq. ([65)) may be written as 

TdS = d[{p + p)V] - Vdp (66) 

The integrability condition is necessary to define a perefect fluid as a thermo- 
dynamical syetem |52]— [55j . It is given by 

dp=P±^dT (67) 

Plugging eq. ([57)) in eq. ([M)) . we have the differential equation 

1 dT 

dS=-d[{p + p)V]-{p+p)V^ (68) 



we rewrite eq ([68)) as 

dS = d 

where c is constant. 

Hence the entropy is defined as 



{p + p)V , " 
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Let the entropy density be s, so that 



S_ ^ p + p ^ (1 + 7)P 
V ~ T ^ T 



(71) 
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where p — and < 7 < 1. 

If we define the entropy density in terms of temprature then the first law of 
thermodynamics may be written as 



d(pV)+-fpdV = {l + p)Td 
which on integration yields 

T = 

where Cq is constant of integration. 
From eqs. ((7T|) and (73), we obtain 
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These equation are not valid for 7 — —1. For the Zel'dovich fluid (7 — 1), we 
get 

T = cop= (75) 
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S — p2 

Co 

^ s ^ p3 ^ T 

Thus the entropy density is proportional to the tempreture. we have 
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where To = so = and Sn — — are constant. 
For radiating fluid (7 = i), we get 



T - p3 



(80) 
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s ~ pi ~ (81) 

Thus the entropy density is proportional to cube of tempreture. 

Now the tempreture, entropy density and entropy of radiating universe is given 

by 

T = Too exp (^^j cosh^ x 



(26 + 1) I a I , 2 / f, — ^(4 - h)L^ MLb , . 3 | a 

s = sqo exp (^-^--^ cosh^^ (Vl " 1^) x 
(25 + 1) I a I , 2 / /] — h{4 - b)L^ MLb , . 3 | a 

S = Soo exp {^0^) cosh^ ( v^t) x 
(?^±1^ tanh^ (^tU ^^4^ + ^ - - ^ 



4(6-2)2 (6-2) 



(82) 



(83) 



(84) 



where Too = soo = ^^^Ve ^^'^ ~ constant. 

5. Case(iii): ab > 

In this case we obtain 

^ = "^^K(^)' ^''^ 
H = I3cosi{y/at + to), (86) 

(i)-2) 

1/ = rcos^^(Vai + to), (87) 
g = exp ( „ + Mx ) , (88) 



Therefore, we have 



2(6-2) 

'=^^^P(2(S)-^^)- ^''^ 



^ = exp _ ^ J cos^^ (x/at + to), (90) 
B = Gexp(^^^^— ^+Mxjcos5(VSt + fo), (91) 



14 



C^Hexp(^-^^!^~Mx^cosi{Vat + to). (92) 
Here E, G and H are already defined in Section 4. 

After using suitable transformation of coordinates, the metric ([Ij reduces to the 
form 

ds^ - exp ( ■^^) cos^(V^r)(dX2 - dT^)+ 



ib-2) 

exp (jy^ + 2MX^ cos{VaT)dY^+ 
exp ( f - 2AfX ) cos(V^r)dZ2, (93) 



,(&-2) 

For the specification of displacement vector f3(t) within the framework of Lyra 
geometry and for realistic models of physical importance, we consider the fol- 
lowing two cases by taking /3 as constant and also /3 as function of time. 

5.1. When 13 is a constant i.e. (5 = [3^ (constant) 

Using Eqs. ([90]) . (|9T|) and (|92l) in Eqs. (fTT|) and ([M)) the expressions for pressure 
p and density p for the model (|93p are given by 



\PI (94) 



1 f 2LX\ 2^ , ^ , 



((3b -A) 2 / ^ N A b(b + 4)L'^ _2 ^'^^^ 



1 f2LX\ ^ , ^ , 



2-b ^ 

' X 



a(3b - 4) 2 / ^.^N ^(4 - 36)2.2 ^^^^ 



1^1 (95) 



From Eq. (|17l) the non-vanishing component F12 of the electromagnetic field 
tensor is obtained as 



n 2ML ^ ( ( Lb \ X\ , ^ , 

= G exp I + 2Mj -| cos(v/^T). (96) 

From above equation it is observed that the electromagnetic field tensor in- 
creases with time. 

The reahty conditions (Ellis [47]) 

{i)p + p>0, {ii)p + 3p> 0, 

lead to 
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a(36-4) 



26 



tan^ (VaT) > 



6(6 - 4)^2 2A/i6 
2(6-2)2 + 2-6 



2A/^ - a, 



(97) 



and 



exp 



/'2LX\ 
[2~b) 



a(36-4) 



tan2(y^T) + 3 



26L2 
(6-2) 



4M2 



4ML6 



(6-2) 



The dominant energy conditions (Hawking and Ellis [48] ) 

ii)p~p>0, iii)p + p>0, 



(98) 



lead to 



and 



'-13^ cos V (yST) > exp 



/2LX\ 
l2~6 j 



62^2 



■ tan2 (VaT) > 



(6-2)2 
2MLb 



2M2 - a. 



(99) 

(100) 



26 - 2(6-2)2 2-6 

The conditions (|108p and (|109p impose the restriction on /3q. 

5.2. When /3 zs a function oft 

In this case to find the explicit value of displacement field /3(i), we assume that 
the fluid obeys an equation of state given by eqs. Here we consider three 

cases of physical interest. 

5.2.1. Empty universe 

Let us consider 7 = in (|3H]). In this case p ^ p = 0. Thus, from Eqs. ([TT|) and 
we obtain 



/2LX\ 



^'^^^ = 3^*'''P (,2^j ^^^^ 



2!,2 



L'b 



(6-2)2 



(101) 



Halford |5] has pointed out that the constant vector displacement field in 
Lyra's geometry plays the role of cosmological constant A in the normal general 
relativistic treatment. From Eq. (|10ip . it is observed that the displacement 
vector /3{t) is a decreasing function of time. 

5.2.2. Zeldovich Universe 

Let us consider 7 = 1. Hence, Eq. (|48p gives p = p. 

Using Eqs. (1101), (HH) and ^ in Eqs. ^ and (O, we obtain 



/3'W = ^exp 



l'2LX\ 



2-b — 

COS *> (VaT) 



M2 



L262 
(6-2)2 



(102) 
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From Eq. (|102p . it is observed that displacement vector fi is decreasing function 
of time. The expressions for pressure p and energy density p are given by 



1 f2LX\ 2^ , ^ , 
Sirp = 87rp = — exp ( j—^ I cos " [VaT) 



fc^tan^(V^T) + | 



(103) 



b{4-b)L^ 3 , MLb 
4(6-2)2 2 6-2 

The reality condition (Ellis [47] ) 

(i)p + p>0, (m)/9 + 3p>0, 

lead to 

(36-4)a 2^ _ . 6(6-4)i2 3 MLb a 

5.2.3. Radiating Universe 

Let us consider 7 = ^. Hence, Eq. p8|) reduces to p = 3p. 

In this case using Eqs. (|90|) - ([92]) in Eqs. (fTTj) - (fT4|) . the expressions for 
I3{t), p and p are obtained as 



2 / 2LX\ 2-b 

■ exp [ - — - ) cos I" (VQ!^^) 



3£;2 " V2- V 



a(36- 2) 
26 



a(56-4) 



46 



tan^ (v^T) 



362L2 



4(6-2)2 



Sirp - 



1 f2LX\ ^ , ^ ^ 
exp I ^— ^ I cos f (VaT^) 



£2 



ML6 
6-2' 



62.2(4 _3&) a . 1 2 ^ ^ . 

^(^-26-2^^ -2"*^^ (^^) 



3 / 2LX , ^ , ^ , 
BTrp = — exp [ y ) cos ' (VQ^^) 



£2 ^^^^ V 2 - 6 



A/L6 
6-2" 



62.2(4 -36) a 
4(6-2)2 - ^ 



1 



2M^ - -atanN VaT 



(105) 



(106) 



(107) 



From Eq. (|105p . it is observed that displacement vector /3 is decreasing function 
of time. The reality conditions (Ellis |47] ) 



(i)/9 + p>0, (m)p + 3p>0, 
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and the dominant energy conditions (Hawking and |48) ) 



{i)p-p>0, {n)p + p>0, 

lead to 

a _ , MLb bL^(i-3b) a o 

-tan^(V^T)>^ + ^A__i___2M^ 

and 

a n. _ . MLb bL^{4:-3b) a o 

- tan^ iV^T) > ^ + _ _ _ 2M^ (108) 

respectively. 

The expressions for the expansion 0, shear scalar cr^, deceleration parameter q 
and proper volume V for the model (|103p are given by 

, _ (2 - 35) exp ( ^) cos^ (V^T) tan (V^T) , (109) 



2bE ^ V 2 - 5 



where C2 = c-^^^'' . 



(J^ = exp (^Y~b) (^/^'^) tai^^ (VaT), (110) 



362 £;2 

86(6 - 1)^2 (i _ tan2(^^T) 
9(36 - 2)2 exp (^) cos2 (VM?') tan^ (v4^r) 

= = i^^ exp (^^1^) cos^ iV^T) (112) 



(111) 



From Eqs. p09p and (jllOp we obtain 

^2 4 



constant. (113) 



612 3(36- 2)2 

The rotation uj is identically zero. 
The Hubble parameter H is given by 

^ = ^-^^§^ exp cos^ iV^T) tan (v/^T), (114) 

For I > 3, the model (j93p starts expanding at T = and attains its maximum 
value at T = ■ After that 9 decreases to attain its minimum negative value at 

T — -r^- Thus, the model oscillates with the period 7ri=- Since § = constant, 
the model does not approach isotropy. The sign of q indicates whether the model 
inflates or not. A positive sign of q corressponds to standard decelerating model 
where as the negative sign — 1 < g < 0, indicates inflations. Recent observation 
shows that the deceleration parameter of universe is in the range — 1 < g < 
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Time (T) 



Figure 2: The plot of deceleration parameter (q) vs. time (T). 



and the present miiverse is undergoing an accelerated expansion [49l [50] . This 
behaviour is clearly shown in Fig. 2. Also the current observation of SNe la 
and CMBR favour an accerating model {q < 0). From Eq. it can be seen 

that the deceleration parameter q < when 



T < 



1 



tan 



(6-2) 



It follows that our models of the universe are consistent with recent observations. 
5.3. Thermodynamical behaviour and entropy of universe 

From eqs. (|75|) . (76) and (71), the expression for tempreture, entropy density 
and entropy for zel'dovich fluid (7 — 1) is given by 



T — Tq exp 



LX 

2-6 



2-b 

cos 2b 



a \T 



b{A-b)L^ 3 
4(6-2)2 2 



(36-4)a 
46 

MLb 



tan^ (VaT) 



a 
2 



So exp 



LX 



2-b 

COS 2b 



(/I 



a IT 



6-2 

(36-4)a 
46 



(115) 



tan^ (V^T) + - 
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S = Sq exp 



6(4 - b)L^ _ 3 
^ 4(6-2)2 - 2 

{b+l)LX\ 3b-2 



2-6 



j cosT^ (^^r) 



6(4 - b)L^ 

4(6-2)2 2"' '6-2 



MLb 
6-2 

(36-4)a 
46 

MLb 



(116) 



tan^ (v^T) + - 



(117) 



where Tq, sq and 5*0 are already defined in section 4.3. 

From eas. (|80|) . (81) and (71), the expression for tempreture, entropy density 
and entropy for radiating fluid (7 = ^) is given by 



T = Too exp 



LX 

4-26 



2-b 

COS 4b 



b(4~b)L^ 3 



(36-4)a; 2 / ^rpN , a 
tan {V^T} + - 



MLb 



4(6 - 2)2 2 6 -2 



(118) 



s = soo exp 



iLX 
4-26 



6-b 

COS 4b 



(36- 4)a 
46 

MLb 



tan^ (^/a^) + - 



6(4-6)£2 3 2 

4(6-2)2 2 ' ^6-2 



(119) 



S — Sqo exp 



/(26+ 1)LX 
V 2(6-2) 



(^^-')"tan^ (v/^r) + | 



6(4-6)L2 



46 

3 



- -Af^ 



AfL6 



4(6-2)2 2"" '6-2 
where Tqq, Sqq and Sqq are already defined in section 4.3. 

6. Case(ii): ab = 

In this case we obtain 

f ( 
f = n exp 

V(&-2) 

/i = {cit + to)K 



(120) 



(I1-2) 
(cii + to) 26 



g — N exp 



Lbx 
2(6-2) 



Afa; 



(121) 
(122) 

(123) 
(124) 
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Therefore, we have 

A = £;oexp(^^^^^(cii + io)^, (126) 

B^Nexp(^^^^+Mxyc,t + to)K (127) 
C = i^oexp(^^^^^-Ma;^(cl^ + to)^ (128) 

(2-h) , 

where Eq = nm b ^ Hq = j^. 

After using suitable transformation of coordinates, the metric ([T|) reduces to the 
form 

ds' = E'„ exp {^^) (ciT) idX' - dT')+ 

exp (^^^ + 2AfX^ (cir)dy2 + exp (^^^ - 2MX^ {ciT)dZ^ (129) 

It is observed that model (129) is same as obtained by Pradhan and Shyam 
Sunder [34]. Thus the physical and geometrical properties of model are similar 
to the model obtained by Pradhan and Shyam sundar [34] . 

7. Discussion and concluding remarks 

In the present study, I have investigated the plane symmetric inhomogeneous 
cosmological models of perfect fluid distribution with electromagnetic field based 
on Lyra's geometry. The source of the magnetic field is due to an electric cur- 
rent produced along the z-axis. The free gravitational field is assumed to be of 
Petrov-type II non degenerate. It is observed that the gauge function I3{t) is 
large in begining and reduces fast with the evolution of universe for all cases. 
Also it is found that f3{t) decreases as time increases therefore I3{t) is decreasing 
function of time and it play the same role in Lyra's geometry as as cosmological 
constant A{t) in general relativity. It means that the displacement vector (3{t) 
coincides with the nature of cosmological constant A. The nontrivial role of 
vacuum in the early universe generates A-term that lead to inflationary phase. 
Therefore the study of cosmological modelsin Lyra's geometry may be rela- 
vant to inflationary models. There seems a good possibility of Lyra's geometry 
to provide a theorectical foundation of relativistic gravitation and cosmology. 
However, the astrophysical bodies is still an open question. In fact, it needs a 
fair trail for experiment. It is seen that solution obtained by Pradhan, Yadav 
and Singh |34| are particular case of my solution. Generally models represent 
expanding, shearing and Petrov type II non degenerate universe in which Sow 
vector is geodetic. Also f 7^ 0, thus models do not approach to isotropy. The 
value of deceleration parameter (q) is found to negative (— 1 < g < 0), imply- 
ing that our universe is accelerating which is supported by SNe la and CMBR 
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observations. This behavior is clearly depicated in Figure 1 and 2. From Fig. 
1, we note that at early stage of universe deceleration parameter (q) oscillates 
between g < and q > —1 and afterwards it will be uniform (negative) for ever. 
This has physical meaning. 

The idea of premordinal magnetism is appealing because it can potentially 
explain all large scale fields seen in the universe todays, especially those found 
in remote proto galaxies. As a result, the literature contains many studies that 
examine the role and implications of magnetic field in cosmology. Maarteens 
[55] in his study explained that magnetic fields are observed not only in stars but 
also in galaxies. In princple, these fields could play a significant role in struc- 
ture formation but also affect the anisotropics in cosmic microwave background 
radiation [CMB]. Since the electric and magnetic fields are interrelated, their 
independent nature disappears when we consider them as time dependance. 
Hence, it would be proper to look upon these fields as a single field - electro- 
magnetic field. It is worth mentioning here that magnetic field affects all the 
physical and kinematical quantities but it does not affect the rate of expansion. 
Also we see that in absence of magnetic field, inhomogeneity of universe dies 
out. This signifies the role of magnetic field. The present study also extend the 
work of Yadav and Bagora |57j with in the framework of Lyra's geometry and 
clarify thermodynamics of plane symmetric universe by introducing the inte- 
grability condition and tempreture. A new general equation of state describing 
the Zel'dovich fiuid and radiating fluid models as a function of tempreture and 
volume is found. The basic equations of thermodynamics for plane symmet- 
ric universe has been deduced which may be useful for better understanding of 
evolution of universe. 
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